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Abstract 

We develop a field theory approach to light propagation in a gas of resonant atoms taking into 
account vector character of light and atom-atom interactions. Within this approach, we calcu- 
late the propagator of the electric field for both short and long-range density-density correlation 
functions of the gas. 
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I. INTRODUCTION 



Transport of waves through strongly scattering disordered media has received much at- 
tention during past years. The character of this process is determined by the multiply elastic 
scattering of light. In a sample of randomly distributed scatterers, the initial direction of 
the wave is fully randomized by multiple scattering, and a diffusion picture seems to be an 
appropriate description of light propagation in the case of small densities of scatters 
Accordance to this picture, he transmission coefficient of a medium is inversely proportional 
to the sample thickness, i.e. follows the familiar Ohm's law. Despite successful predictions, 
this theory does not take into account interference effects. P.W. Anderson was the first 
to demonstrate that sufficiently strong disorder can result in a localization of quantum- 
mechanical wave functions [2]. This phenomenon is referred to as Anderson localization. In 
fact, even away from the localized phase, in the so-called weak localization regime, wave in- 
terference affects the physics of light propagation. For example, the constructive interference 
of counter-propagating light amplitudes results in the coherent backscattering peak 0, J]. 
In the past 20 years, mesoscopic physics has been developed, with beautiful experimental 
and theoretical results pertaining to the weak localization regime 5), (], Q|. 

A system of randomly distributed atoms with an atomic transition frequency close to the 
light frequency, provides necessary conditions for experimental observation of backscatter- 
ing |s, U 10, In fact, under these conditions, light scattering on an individual atom 
is characterized by a large cross-section. However, in contrast to classical scatters, in this 
system one has to take into account resonance dipole-dipole interactions that alter the char- 
acter of the light propagation. There are two different mechanisms of excitation transfer. 
According to the first mechanism, an excited atom decays into the ground state by emitting 
a photon. Then the photon propagating through a system of resonance scatters, is absorbed 
by another atom, and thus excitation transfer takes place. In the second mechanism, an 
excitation may transfer from atom to atom nonradiatively by the resonance dipole-dipole 
interaction. This is equivalent to including the longitudinal component of the electromag- 
netic field. This mechanism becomes more important with increasing the density of atoms 
n. If the average number of atoms per cube with the size of the wavelength A exceeds unity, 
n\ 3 > 1, the typical energy of dipole-dipole interactions is comparable with the inverse of 
the the radiative lifetime of excited atoms, and the transport of nonradiative photons dom- 
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mates the transport of radiative ones. The resonance dipole-dipole interaction a 



important role in the radiative transfer by binary collisions of atoms in gases 
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so plays an 



id- 
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This problem has been actively studied under in the case of low atom number densities 10|, 



15] . By using the scalar model of light it has been demonstrated that interference effects 



are substantially suppressed as a result of the motion of scatters [16|]. In Refs. [ll| the authors 
considered the vector character of light and took into account atom-atom interaction. It was 
found that close to the resonance, the transport slows down considerably because the light 
is captured by atoms for a longer time. Such a behavior results in a considerable reduction 
of the speed of light. This effect have been observed for the classical systems, but for the 
dipolar gas the effect differs [§, [9]. 

In this paper we focus our attention on the case n\ 3 <C 1 and consider different density- 
density correlation functions for light scatters. The latter will be treated as a (1 + 3)-level 
atom (In other words, we consider atoms with a resonant s—p transition.). We also take into 
account the effects of dipolar atom-atom interactions on the light propagation. The influence 
of the atomic density-density correlation function ( the second cumulant, or second central 
moment of the probability distribution for different atomic realizations) is of special interest 
for us. It appears that in the considered case of a low density system, only this correlation 
function affects the distribution of total transmission coefficient. We therefore assume that 
the other cumulants are zero (that is we consider gaussian probability distributions). This 
means that a general high order density correlation function reduces to the product of pair 
density-density correlations. 

The paper is organized as follows. In Sec. II, we introduce the action for the (1 + 3)- 
level dipolar gas interacting with quantum electromagnetic field. An effective action for 
the light propagating in the gas, is derived in Sec. III. We evaluate the propagator for the 
electric field and identify the localization length of light in the gas in Sec. IV. There we also 
present simple relationship between the density-density correlation function and the light 
propagator, and analyze the result for various correlations in the gas. Finally, in Sec. V we 
conclude. 

Throughout the paper we adopt the system of units in which the speed of light c, the 
electron charge e, and the Planck constant h are set to unity, c = e = h = 1. 
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II. FIELD THEORY DESCRIPTION OF THE SYSTEM 



A. General considerations 

The propagation of light in a medium is described by the Maxwell equations supplemented 
with ones characterizing the properties of the medium. We express the solution in terms 
of propagators that relate the amplitude of the electromagnetic 4- vector potential A k (X) 
(k = 0, 1, 2, 3) at a space-time position X = t, x, with A n (X') and is defined as 

D kn (X - X') = (G\TA k (X)A n {X')\G) (1) 

where \G) is the ground state of the system, when the all atoms in the ground states and 
the photons are absent. We have used above the fact that the medium is transactionally 
invariant. For the medium that consists of scatters at random positions, the propagator of 
a given sample depends on the location of scatterers. To get general characteristics for light 
propagation in the disorder medium one have to average the propagator over disorder: 

D kn {X-X') = ((G\TA k (X)A n (X')\G)) dls (2) 

The averaged propagator is referred to as the amplitude propagators. 

For further calculation the quantum field theory methods are applied. Thus, for the 
evaluation the expression ([2]) we have to construct the action for the light interacting with 
dipole gas. The action for the system depends on atoms degrees of freedom as well as 
electromagnetic field one. In a general form it can be written as 

S = S(A k ,^,W) (3) 

where are fields corresponds to atoms and A k 4- vector of electromagnetic field. Therefore 
the partition function Z for the system is 

Z=j J] VA k Y[v{W,W}exp[iS(A k ,WJ?)] (4) 

"* fc=0,l,2,3 

On the base of the participation function (jlj) we can evaluate the properties of a certain 
realization of disorder. To find the averaged propagators (J2j) one have averaged over disorder. 
Applying continual approach we figure distribution for realization of the number density p(x) 
as: 

P[p(x,t)\ = Cexp 



n „ 

/ dx\ . . . dx. 



p{x 1 ti)p(x 2 t 2 ) 



n!r( n )(xi 



x n , t\ 



t n ) 



(5) 
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Here V( n >(xi . . . x n ,ti . . .t n ) is n-cumulant for the distribution. The first cumulant is equal 
to the averaged number density p (x,t), while the second one coincides with dispersion for 
the number density : 

r (2) (x]ti, x 2 t 2 ) = {{p{xib x ) - p (x 1 t 1 ))(p(x 2 t 2 ) - Po{x 2 t 2 ))} (6) 

If the number density of scatters is not high, the other cumulant are close to zero. For 
further evaluation we take into account that the number density can evolve with time, then 
the improved distribution ((Sj) for p(x, t) can be depicted as: 

P[p(x, t)] ~ exp 

There the averaging-out over disorder realization results in appearing a supplementary field 

p(x,t) 

Z= f J] VA k ]Jv{W,$}Vp(x,t)P[p(x,t)]exp[iS(A k ,W,^ (8) 

fc=0,l,2,3 

B. Model for atoms. 

For the evaluation the field propagator we are in need of the full action for the system 
S(Ak,^,^2)- We start the analysis from consideration of a single atom immersed into an 
external electromagnetic field. The total action for a single atom S to t_atom in the electromag- 
netic field consists of three terms: the action of the free atom S a t om , the free electromagnetic 
field Sf and the term caused by the interaction of the electromagnetic field with the atom 

Sint- 

Let us find the action of a free atom. All evaluation will be carried out in the assumptions, 
that the atom is infinitively heavy. It allows us to take into account the evolution of just 
internal degrees of freedom. The translation degrees of freedom of the atom are not taken 
into account, because the effect of the momentum exchange between the atom and incident 
electromagnetic field is supposed to be small. In general case atomic structure is fairly 
complex and with increasing the number of electrons in the atom its complexity is increasing 
and the number of energy states grows [l?], ■ The action of the external electromagnetic 
field inspires transitions between the atomic states. If the atom with ground state energy 
e g and excited state energy e a is subjected to the monochromatic radiation with frequency 



,(p(x,t) - p (x,t)) (p(x',t r ) - pojxXj) 
T( 2 \x,t;x',t ; ) 



(7) 
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uj and the conditions uj ~ |e — e g | holds true, the most of transitions occurs between the 
ground and the excited state. The others states does not affect on the atom dynamics in 
the electromagnetic field. It is assumed the character of the resonance transitions to dipole. 
In accordance with the selection rules for the dipole transitions, the final orbital momentum 
for an atom Lf and initial orbital momentum Li are related by |Lj — Lf \ = 1. Thus, if the 
ground state is singlet, the excited state has to be triplet, then transitions takes between 
four states \E, I) with the energy E and projection of orbital momentum I 

|^> = |e fl ,Z = 0> 
R = k.,i = -1> (9) 

|0) = |£ a ,/ = 0) 

|+) = \e a ,l = +1) 

So consider the system that contains atoms in the lowest (E — e g , I — 0) and first excited 
states (E — e a , I — —1,0,1) where we let cf with % = g,—,0,+: create these respective 
states. If \ Vac) is an atom vacuum state (i.e. the state that corresponds to the atom absence), 
these operators satisfy the following condition: 

\i) = cj\Vac) \Vac) = Ci\i) i = -1,0,1, g 
{c+, c k ) = 5 ik {c+, cl) = 0, {a, c k } = 0, Vi, k 

We will use further the change of basis c+,c+ = (cl + c+) /y/2 , c y = 
— i (c+ — cl) /a/2 , c+ = Cq, and the hamiltonian of a free atom can be described in 
new terms of annihilation and creation operators q, cf where i — g,x,y, z as 

H £aC x ^aCy Cy -\- S a C z C z -\- £gC g Cg (H) 

Let us introduce the Grassmann fields ipi and ^ that correspond to annihilation and cre- 
ation operators Q,c+, i = g,x,y,z. If the fields ipi and ^ correspond to coordinate and 
momentum, then the free atom Lagrangian takes the form: 

L a tom= $i(t) id Mt) - e g 4> g (t)4> g (t) - E a ^ i>a(t)lpa(t) (12) 

i=g,x,y,z a=x,y,z 

and the action for a single atom 



'atom 



-I 



dt 



$i(t)idti/>i(t) - e g i) g {t)^ g {t) - e a ^ 4> a (t)i> a (t)+i5 £ A(t)^(t) 

i=g,x,y,z a=x,y,z i=g,x,y,z 

(13) 



A short imaginary parts appears to satisfy causality principle for atom evolution. Taking the 
energy of the ground state is equal to zero we can rewrite action for a single atom as 



Satom = J ^g(t) [ify + iS\ ^ g (t)dt + J $ a (t) [id t - A + l5\ ^) a {t)dt, A = S a - S g (14) 

C. Interaction of atoms with light 

The other two terms of the total action (S to t_atom) Sf and S int depend on the properties 
of electromagnetic field. If the electromagnetic field is determined by 4-potential Ai, the 
action associated with the free e-m field fl9l. 0] is 

S f = ~— I d A x ( -F ik F ik + -(d k A k ) 2 ) F ik = ^ - ^ (15) 
i 4ttc J V 4 2C k! J lk dx l dx k K ' 



The constant £ appears as result of gauge invariance for the electromagnetic fields |19l. |20 |. 
Further we choose an appropriate gauge for evaluations. The term resulted from the inter- 
action has to consists of field variable Ai as well as 4-current j l , which is determined by 
atom properties: 

S int = ~^J d A xA l3 l (16) 

If the atom is placed at the point x n and the wavelength of electromagnetic field surpasses 
remarkably its size, we can treat the atom as point dipole and the 4-current produced by it 
is 

n xt) = { cp } xt) 1 = ( 1 (17) 

[ j(xt) j I di(t)8(x — x a ) J 
Integrating expression ( fl6l) over time and applying for the 4-current (fTTl) . we get : 

S int = J dAxA if = \J E *( x n, t)d a (t)dt a = 1, 2, 3 (18) 

where E a (x,t) is electric field strength: 

E a (x, t) = - [d a A (x, t) + d A a (x, t)} (19) 

If do is the value of dipole moment associated with transition e g — > e a , the component of 
the operator d a (t) is written as 

d a {t) = d [^ a (t)^ g (t) + $ g (t)lf> a {t)] 
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The expression ( (TBI) is exact for point-like dipole. 

Thus the total action for the system of the atom and the electromagnetic field: 

Stotatom = ~^- c J d^x {^\F ik F lk + ■^(d k A k ) 2 ^ 

\ J dtE a (x a , t)do [Tp a (t)ip g (t) + $ g (t)il> a (tj\ (20) 
/ dt$ g {t) [id t + i5] il> g (t)dt + / $ a {t) [idt — A + iS] i[) a (t) 
The total action of the atoms ensemble Stotai and the electromagnetic field can be found by 
the generalizing (1201) . Naively one may assume the total action for an atomic ensemble is 
sums actions of free atoms and interaction terms over all ensemble, and the free electromag- 
netic field term. Applying expression for S% tom the action of a single atom at a point x n (TTij) 
and its current J l (x n ) (ITT)) : 

S' M al = £ S »tom -- c f A *H ^ + S f W ( 21 ) 

n n 

But this naive approach leaves out of account the interactions between atoms, resulted in 
an additional term in the total action. The force acted between two atoms is referred to as 
the Van der Waals force. The elementary processes on which the Van der Waals interaction 
is based is the exchange of the pair photons between the two atoms. This force depends on 
the distance as well as energy of resonance transitions. If the distance between the dipoles 
R is short, the virtual photons propagate instantly and the force is proportional to R~ 7 . In 
the opposite case, the delay in propagation of photons should be taken into account. The 
character dependence force on the distance is changed and has a form R~ 8 . The interaction 
between dipoles has electromagnetic nature. Thus the electromagnetic filed in the dipole 
gas consists of two parts: the first one originates from the external source A ext and the 
second part appears as a result of interaction between dipoles Af 1 . Further we are refereed 
to A^ xt ,A^ n as "external" and "quantum" fields. Therefore the total action for dipole gas 
in the electromagnetic field takes the form: 

Stotai = E S atom ~ \ J (A^ + A T) £ '"M + S, {Af) + Sf {A?) (22) 
n n 

If the external field is equal to zero, one can integrates over the quantum field A q and the 
new effective action is determined by the atomic freedom of degree. This action includes sum 
free atomic action for all atoms [13] and the term Saa related with interquartile interaction 

S AA = -\ [ d^x' Jm ^ Jn ^f ) 8{r -r' - \x-tf\) (23) 
c 6 J x — x'\ 



S 



where J m = Xli m ( x n) is the total 4-current caused by all dipoles. If we consider the limit 

n 

case, when the interaction can be treated as instant, i.e. A • R mn jc <C 1 [R mn is distance 
between atoms n and m) the interaction takes the form: 



SaA — — — j ' Hrnn{t)dt 



m,n,m^n 



da (t) d™ (t) R mn — 3d™ (t) d m (t) R a ,mnR/3,mn 

R~ 5 

i.e. we get the dipole-dipole interaction for static limit, as we should. The total action S tota i 
for the dipole gas is given by surri22l with i*(x n ) determined by [T71 

Because the ground state of the system can be determined by as product of the atomic 
the ensemble ground state and the electric field ground state |0) (the state when the number 
of photons is equal to 0), the ground of the system is 

iG>=io>ni0>»=i°> iri fqc >« 

n n 

the partition function for the ensemble is 

Z= lim lim Af [ T>\A extfi A ext > 1 A ext ' 2 A ext ' 3 }V\A qn '°A qn ' 1 A qn ' 2 A qn ' 3 } 

«i->-aot/-«x> J 1 ' 1 ' 

\{ \{ V^V^(t f )^(t t )-eMiS total ] 

n «=g,l,2,3 

here TV is fields independent normalization constant. 



III. EFFECTIVE FIELD THEORY FOR LIGHT 



Above we have found the partition function for the interacting atomic ensemble. The 
paper object is to derive the Greens function for electromagnetic field. To carry out this 
deducing, we have to integrate over atomic fields ip™. The integration over atomic variables 
in the partitions function for the interacting dipole gas can not be performed accurately and 
some approximations must be employed. For weak interaction, perturbation theory can be 
applied, when the interacting term exp (iSi n t s ) is expanded in the form of a power series: 

Cn2 qn4 

exp(zS ints ) = l-^ + ^ + ... 
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we keep the terms up to the forth order. For brevity it is convenient to introduce the 
following notation for the fields: 

(A qn Y + (A ext Y 



A 1 = 



B l 



V2 

(A" n y - (A ext y 
7! 



Then the partition function for the dipole gas is 



Z — Z\ ■ Z 2 



n [v^v^}^(t f )^(t i )(i-^+^) 

i=g,l,2,3 \ - J 

Z x =M A \v { A A 1 A 2 A 3 } exp[^ e J J] Z n 
Z 2 =N B jv {B°B 1 B 2 B 3 } exp[i^ e J 



where 



S fieid - 4^f d4x 



02 

Afield 4?rc 
V2dp 



\F m F™ + ±{d k A k Y 
±P tm P™ + 7 L(d k B k y 



P _ 9Ak _ OA 



dx l 



dx k 



1 tk dx* 



dB k dB % 



dx k 



int 



J dr ■ E a (xr) [C(r)^(r) + ^(r)^(r)] 



Satom — I dT 



exp{iS™ 



atom} 



(24) 



E CM (idr -A + iS) C(r) + r g {r) (id T - iS) ^(r) 

a=l 



E a (x n , r) = - (d a A (x n , r) + d A a (x n , r)) 
and NaiNbiNu are constants. Because the atomic fields do not mix in full action, the 
integration for each atom can be performed independently. As result of integration we get 
finally for Z n 

( l-i 2 -^J dt x dt 2 E a {x n , t x )E a (x n , y • G e (h - t 2 ) N 

-i i^f) 2 / dhdhdhdU ■ Ge (*i - h) G g (t 2 - t 3 ) G e (t 3 - U) (25) 
y ■E a (x n ,t 1 )E a (x n ,t 2 )E h (x n ,t z )E h (x n ,t A ) ) 

where the atomic Greens function for the ground G g (t 2 — ti) and the excited states 
G e (h — h) are defined as: 

(^(t 2 )^ n (ti)) = iG g (t 2 -h) = e (t 2 - 1,) 

<C(*2)^?(*i)) = iSabGe (t 2 -U)=e {t 2 - h) exp [-iA (t 2 - h)] 5 ab 
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The final expression for Y[ Z n up to the forth order for do is 

n 



( 1 - Yj^ I dt i dt 2Ea{x n , t x )E a {x n , t 2 ) ■ G e (h - t 2 ) ^ 

n 

+1 Y 5 dt x dt 2 dt 3 dt r G e (t x - t 2 ) G e (t 3 - t 4 ) 

-E'a(-^'n) t\)Ej a \X n: t^Ebi^Xm, t^E b (^X m , t^j 

~Y l {if)* J dtidt 2 dt 3 dU ■ G e (h - 1 2 ) G g (t 2 - t 3 ) G e (t 3 - t A ) 



(26) 



i -£< a (x n , t-^)£' a (a; n , t 2 ^)E b {x n , t 3 )E b yX nj t^j 

If we introduce the atom number density as 

P( x ) = Y S ( X ~ Xn " > 

n 

The discrete description of atomic ensemble can be replaced by continuous one and equation 
transforms into 



/ 



dhdt 2 d 3 x ■ p (x) E a (x, t x )E a (x, t 2 ) ■ G e (ti - t 2 ) 
i-^j f dtidt 2 dt 3 dt i -d 3 xd 3 y- [p (x)p (y) — p (x) 5 3 {x — y)~\ 
■G e {h - t 2 ) G e (t 3 - i 4 ) E a (x, t 1 )E a {x, t 2 )E b (y, t 3 )E b (y, t 4 ) 
/ | f dt\dt 2 dt 3 dtt±d 3 x ■ p[x) ■ G e (ti — t 2 ) G g (t 2 — t 3 ) G e (t 3 — t 4 ) 
■E a (x, t l )E a (x, t 2 )E b (x, t 3 )E b (x, t 4 ) 



\ 



\ 



J 

(27) 

The dipole density p(x) is stochastic value and depends on realization of disorder in 
the investigated sample. Unless the problem under consideration is light propagation in a 
certain sample, the answer has to averaged over density realization. Denote by P [p{x)} the 
probability density that the atom number density is equal to p(x). The rigorous definition 
of this density reads 

P[p(x)] = Cex V [L [p(x)}} 
where the normalization constant is defined as 



C- l = J D P (x)exp[L [p(x)}} 

Generally, realization of disorder may vary with time (it means, in the course of evolution 
the corresponding cumulants change their values), but we restrict our investigation by the 
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frozen disorder, where the first and second cumulant are not equal to zero. This approxima- 
tion holds true for low density, when a typical contributions of the higher cumulants drops 
remarkably with increasing its order 2l| . so L \p(x)\ can expressed as 



2 J rw{xi,x2) 

Moreover, the atomic gas is assumed to be translation invariant and the density-density 
correlation can be expressed as 

r^(x, y) = (p{x)p{y)) -(% = p 2 C(x - y) + p 5 3 (x - y) 
(p(x)) = p 

where C(x — y) is dimensionless correlation function for density, ( density- density correlation 
is supposed to be independent on the atoms positions and interatomic distance affects on 
the the correlation function ). So averaging over position we get eventually 



\{z n =tiK(t,E,C,p Q ) 



(28) 



/ 



Af 



■2dl 



J dtidt 2 d 3 x ■ E a (x, t 1 )E a (x, t 2 ) ■ G e (ti - t 2 ) 



\ 



A 2d 0P0 



J dtidt2dt^dt4 1 -d 3 xd 3 y 
■G e (h - t 2 ) G e (t 3 - U) E a (x, t x )E a {x, t 2 )E b (y, t 3 )E b (y, * 4 ) 
+ 1 2 j dt 1 dt 2 dt 3 dU-d 3 xd 3 yC{x - y) 



■G e - t 2 ) G e (t 3 - h) E a {x, tjEaix, t 2 )E b (y, t 3 )E b {y, t 
—i f dtidt 2 dt 3 dt 4: d 3 x ■ p ■ G e (t 1 — t 2 ) G g {t 2 — t 3 ) G e (t 

■E a (x, t 1 )E a (x, t 2 )E b (x, t 3 )E b (x, t A ) 



-U 



(29) 



and the partition function is given by 



Z X =M J V {A°A l A 2 A 3 }exv[iS} ield + \ogK.{t,E,C,p Q )\ 

IV. PROPAGATOR FOR THE ELECTRIC FIELD 

In the previous section we have deduced the partition function for the atomic ensem- 
ble. On the basis of it we can establish a character of evolution for 4-th potential A^ (xt) , 
characterized the electromagnetic field evolution. In itself the 4-th potential is not physical 
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observable value, so we focus our attention on electric field and calculate the propagators 
for the electric field: 

Da^RihRzh) = -% {E a {R x t x )E b {R 2 t 2 )) 

The propagators of Greens function carry out complete information about spectrum or dy- 
namics of electric fields. It describes the propagation of waves form an external 5-delta 
function source. For further derivation it is convenient to express the action in terms 
of the electric field. Namely, the 4-vector A a can be expressed in terms of 4-vector in- 
cluded components of electric field strength E a (xr) supplemented with zero component 
E (xr) = — d A — d a A a . Hence, the action of free electromagnetic field!5l in impulse- 
frequency representation converts into 



d 3 k dtt 



-E (k, Q)E (-k, -SI) + E a {kQ)E b {-k, -Q)U ab {k, Q) 



S 11 = — 

field g7rc j 2n l £ 

tt n <~>\ kakb r ^ 2 k\ 
LLab{k, il) = -Jp- + ab 

The zero component of 4-vector (Eq, E a ) ,a = x,y, z corresponds to the calibration and does 
not influence on any physical observable values. This fact is proved by explicit representation 
of Sj ield , and fC (t, E, C, po) does not contain E . 

The Greens function of the photon state is then determined by 

#2*2) = -i {E a {R x ,T x )E h {R 2 T 2 )) 

AT / n VE k exp [iS) ield (E k )} K (E k , t) ■ E a (R v TJE^R^) 

fc=0,l,2,3 

#J EI VE k exp[tS} ield (E k )])C(E k ,t) 

fe=0,l,2,3 

this formula can be represented diagrammatically(see fig[T]). 

The right hand side features the contributions of free evolution (diagram A), single scat- 
tering (diagram B), two scattering on two atoms (diagram A) and double scattering on the 
same atom (diagrams D-G). The explicit expression of the free photon propagator is 



D2 6 (fl 1 r li i2 2 T 2 )= -i (E^T^E^R^)) 

f d 3 kdfl Aixc 
= J 4 exp [-iO, (T 2 - T{) +ik c (R 2c - Ri c )] ^ 2 - fe2 - ^ • (fl 2 5 ab - k a k b ) 
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FIG. 1: Diagram representation for the photon Greens function. The wavy, solid and dash lines 
correspond to free photon propagator, the ground and excited atom Green functions respectively. 

while the first order contribution is given by: 

Dl^RiT.R^) = J dt x dt 2 d?x- [G e (h - t 2 ) + G e (t 2 - h)} 

■ D^R^xtjD^xt^n) 

The second order term contains two different contributions. The first one appears as the 
second order Born term for the one-atom scattering and can be presented as: 



D^iR^RiT,) 



2d 



2\ 2 



p / dtidt 2 dt 3 dt A d?x 



■ D° pa (R 1 -x, Tl - tl ) Dl (x - R 2 , t 2 - T 2 ) ■ D » { °f h) 



j +G e (U - ti) [G g (ti - t 3 ) G e (t 3 - t 2 ) + G g (h - t 2 ) G e (t 2 - t 3 )} \ 
+G e (h - U) [G g (t 4 - t 3 ) G e (t 3 - t 2 ) + G g (t 4 - t 2 ) G e (t 2 - t 3 )} 
+G e (t 4 - t 2 ) [G g (t 2 - t 3 ) G e (t 3 - h) + Gg (t 2 - h) G e (h - t 3 )\ 
V +G e (t 2 - U) [G g (t 4 - t 3 ) G e (t 3 - ti) + Gg (t 4 - ti) G e (tx - t 3 )} j 
The other term is in charge for the consequences photon scattering on two different atoms. So, 
in contrast to the previous the second order term it is proportional to the square of atom 
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density and is equal to 



D 2 J(R 1 T 1 R 2 T 2 ) = (^^j J dt l( lt 2 dt 3 dU-d 3 xd 3 y 



1 ' \ | C(x-y) + C(y-x) 



■ (G e (t 4 - *0 + G e (t, - U)) (G e (t 2 - t 3 ) + G e (t 3 - t 2 )) 

■ D° pq (x -y,U- h) ■ D° ap (R, -x,T x - h) D° bq (y - R 2 , t 2 - T 2 ) 

The straightforward result can be obtained in the momentum-frequency representation, 
where the free propagator is given by 

4ttc 

2 2 i 

The terms D ab and D^ b can be easily calculated and are equal to 

D\ h (k, n) = D° ap (k, n) • x (n) • D° pb (k, n) 

-Daf (*, n) = D° ap (k, n) • x (n) K 9 (fc, n) + (*, n)] x (n) • £>? 6 (*, n) 

4Ad 2 p 1 



x(n) 



c 2 ft 2 - A 2 + i<J 

^ n /,/n\ C {k' - k) + C (k' + k) 



H ab (k,n)= D ab (k',n) 

J (2tt) 



2 1 

The expression for Dj, after the Fourier transformation takes the form 



(*, n) = i [-^ ) poB (n) (*, n) i& (*, n) 

V (fix, ft) = G 9 - Q) (G e (-ft) + G e (fix)) 2 

The integral for i? (fix) is divergent and the singularities are connected to the large-k be- 
havior of the Green's function in Fourier space. So, in order to remove the singularities, 
one can modify this behavior. The standard method, often used in quantum field theory, 
is to introduce a large-momentum cutoff. Because of the vector character of the photon 
propagator, we allow for two regularization procedures with different parameters, and 
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=1 



d 3 k 

d s k 



K (k, SI) 
D° bb (k,n)-Anc(l- 2 -^ 



(2n) 3 

cA| - cSl 2 A T - i2c\n 3 



+ 



d 3 k A ( 2S1 2 

t47TC 1 — 

(2vr) 3 V k 2 



The employed regularization procedure is of course not unique, but it suffices for our pur- 
poses. The physical result is largely insensitive to such details. The final expression for 

d 2 J (K n) 



2?i 1 (fc,n)=i 1 l2 ( f 



°- p B (fl) D° ap (k, n) D° pb (k, SI) 



B (SI) = -id 



2A 



SI 2 - A 2 + i<J 



-id 



2A 



n 2 



fi 2 - A 2 + iS 



2\Sl 3 



— c- 



2A 



SI 2 - A 2 + i(J 



where the constant Ci and C 2 are related with the cutoff parameter. The total propagator 
is given by: 

D ab (k, SI) = D° ab (k, SI) + D° ap (k, SI) ■ X (fl) • -Dp b (fc, SI) (30) 

+ n) • x (n) (*, n) + ^ (*, n)] x (n) • ^ (fc, n) 
+ (*, n) x (n) 



-d^l x (Q)^) .Di(fc,n) 



3p 



4Arf 2 p 



1 



~2 fi 2 _ A 2 + f( J 



/(fib' 
7^-3D ab (k',Sl) 
27T 



C (k' - k) + C (k' + k) 



(31) 



(2vr) d 2 

We have omitted the terms included renormalization constants Ci,d, because they result 
in shifting of the values of dipole momentum do and energy A. 

To specify the problem completely one needs to give the value for the density-density cor- 
relation C (r). We start consideration from white- noise correlation, where C (r) = K5 3 (r). 
Such a function describes short-range correlation appears between the neighboring atoms. 
With accuracy of renormalization constant 



H ab (k,Sl) = -2ic\Sf\K 
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and 



D ab (fc, ft) = D» ab (k, ft) + (k, si) ■ x (n) • ^ (fc, ft) 



+ d° p ft) • x (n) • ^ (*, n) • x (n) • ^ (*, n) 

2d |f 
3p 



— 2r? I ft 3 1 

+ ^ (*, n) x (n) • q (i + ^po) x (n) • D° pb (k, ft) 



x(ftH 4Ac ^ 1 



c 2 ft 2 - A 2 + i<J 

For the low density of scatters with the same accuracy the answer can be present in a 
compact form 



or 



D ab (k, ft) = l>°, (fc, ft) + d% (k, ft) • x (n) • ^ (*, n) • x (n) • D° qb (k, ft) 

2a|ft 3 |(l + fTp )' 



+ (*, n) 



1 



[x(^)]" 1 



3p [x (^r 1 



D° pb (k, ft) 



~ i& (*, ft) + -D° p (*, ft) • x (n) • D° pb (k, ft) 
+ ft) • x (n) • ^ (*, n) • x (n) • ^ (*, n) 

= 4Arf 2 p 1 

c 2 ft 2 - A 2 + f Ad 2 |ft 3 | (1 + Kp ) 



D a6 (ft, fc) = 4vrc 
a (ft) = 



(Uft 2 (l-a(ft))-M; fe 1 
ft 2 (1-a (ft)) - A; 2 + ie (1-a (ft)) 

167rArf 2 p 1 



ft 2 



A 2 + 8 -^^ (1 + K Po ) 



Here a (ft) coincides with permittivity for the dipole gas. Because of the correlation the 
imaginary part of permittivity has increased in contrast to the common dipolar gas. Such 
a behavior leads to increasing the scattering mean free path Ir. It is determined by the 
imaginary part of the pole 



The result reads 



l = Im^ (1 _ ft(n)) 
Mr 



2 |ft 2 | (1 + K Po ) (ft 2 - A 2 ) 2 + (r |ft 3 |) 



3tt 



Po 



(r|ft 3 |) z 
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As expected, this result indicates that with increasing correlation scattering is determined 
by multiply scattering, while for the low values K the imaginary part is determined by the 
evolution from the internal degrees. 

Above we have analyzed the short-correlation case, but the situation becomes drastically 
different when long range correlation sets in due to induced particle aggregation, as for 
instance in gelly or aerogels. If it is the case the seried30l can collected as follows: 



D ab (k, n) ~ D° ap (k, n) + D° ap (k, n) % (n) D° pb (k, n) 

+ D° ap (k, n) ■ x (n) [D° pq (k, n) + H pq (k, a)] % (n) ■ ^ (*, a) 



c> w-A*+£&&\m (32) 



AAd 2 oPo 1_ 

3c ^"0 

Thus, with the accuracy upto the second order the electric field propagator can be rewritten 

as 



D ab (k, n) ~ (*, n) + Dg' (fc, n) ■ x (n) H pq (k, n) x (fi) • n) 
D e// fife n) - irr Jg^(i-4ncx(n))-k a k b 1 

The first term of the Green function can be attribute to a photon free propagation in the 
atomic resonance gas without correlation, while the second term is related with "scatter- 
ing on density correlation" in such a medium. Taking into account, that the medium is 
translation invariant, one can decompose H ab (k, Q) as follows 



H ab (*, n) = J ^D ab (k', n) • c(k ' + + = a (*, n) 5 afe + s (*, n) ^ (34) 

Thus, in contrast to the short-range correlation, the long-range density correlations result in 
space anisotropy for the light propagation. The character of this anisotropy is determined 
by the coefficients A (k, Q) , B (k, Q). The function H ab (k, Q) can be rewritten as 

Hob (k, O) = -J cC (R) cos [kR] • {5 ab n 2 + d a d b ) 6XP ^ R] ■ d 3 R 

= -J cC (R) cos [kR] ■ (s ab Q 2 + 6 ab j + d a d b - S ab ^j eXp[ ^ |i?] ■ d*R 
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The tensor H ab (k, Q) can be decomposed into terms: 

H ab (k,Q) = Hl b (k,n) + H 2 ab (k,n) 

Hl b (k, {l) = -JcC (R) cos [kR] • [d a d b - 5 ab ^j eXP ^ |/t>] • d'R 
H 2 ab (k, n) = -JcC (R) cos [kR] • 5 ab (n 2 + |) 6Xp[ ^ |i?] • d'R 
The trace for the first term is equal to 0, so a general expression for H^ b (k, Q) 

H 1 ab (k,n) = A(k,n) (^ ab -3^j 

and the value for A (k, Q) can be found if we multiply the tensor on so 



k a k b i 



k 2 



Ht b (k,n) = -2A(k,n) 



kah f cC(R)co S [kR].(^- 6 -^)R 



k 2 J v ' \ R 2 3 

d ( 1 d ( exp [i \Q\ R] > > . .. rj 
' d R 



dR \RdR \ R 



- j \c (R) cos [kR cos 0] • ^cos 2 
d f 1 d f exp [i \Q\ R] 



dR \RdR \ R 

So after the simplification we get relationship for A (k, fl) 

2sina 
P (a) = — 



2nR z dRd cos 9 



The evolution the second term is simpler, because it is proportional to 5 ab . After the re- 
placement Laplacian we get 

H 2 ab (k, fl) = -5 ab J cC (R) cos [kR] • _ ^(j^ . # R 

Aire 8-jtcQ, 2 f 

= 5 ab —C (0) - S ab —— C{R)- exp [i \Q\ R] sin kR ■ dR 
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The final expression for H a b (k,Q) is 

H ab (k, n) 



Aire 

A(k, tt)-B(k, tt) + C(0) 



6 ab - 3A (k, fl) ^ (35) 



am - -c/ c<«> (™ + • A ( - («)) . ^ 

B (k, n) = f C(R)- exp [i \Q\ R] sin kR ■ dR 

ok J 



From the structure of the expression (j3o]) one can conclude, the only long-range correlations 
cause space anisotropy for photon propagation. 

V. CONCLUSIONS 

We have studied the photon propagation through randomly distributed resonant atoms. 
Starting from the exact action included interatomic dynamics, atom-atom interaction(in 
electric-dipole approximation) and electric field atom interaction, we have derived the Green 
function for the dressed photon that describes both the inelastic processes like a spontaneous 
emission of the atoms and elastic scattering photon in gas with density-density correlations 
via resonant atom-field interaction. The proposed model of atoms allowed us to describe 



consequence of transitions in real atoms and abandoned two-level models 22J mostly used. 
The applied approach are not restricted to the rotating-wave approximation. We have 
shown that the energy associated with transition between resonance states does not change 
for low atomic density where the coupling parameter is 167r ^° Po is small. There are two 
energy-transfer channels in general: one channel through resonant dipole-dipole interaction 
mediated by virtual-photon creation and destruction and the other one through emission 
and absorption of real photons. In particular for strong atom-field interaction, the (over 
a period averaged) energy transfer can be inhibited due to destructive interference of the 
two available transfer channels. The developed approach allows to consider both these 
channels simultaneously without separation. Short-range and long-range correlation for the 
atom densities have to be studied. From the general consideration, it was found that even 
for translation invariant medium, long-range correlation results to space anisotropy for the 
photon propagation. Clearly, the present analysis has left a number of open questions, on 
which future work will concentrate. In particular, the problem of energy transfer in ensemble 
of the moving atoms or influence of saturation effects on the character light propagation. 
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